On Constructing the Minimal Seed of Turbulence: 
Nonlinear Transient Growth 
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Linear transient growth analysis is commonly used to suggest the structure of disturbances which 
are particularly efficient in triggering transition to turbulence in shear flows. We demonstrate that 
the addition of nonlinearity to the analysis can substantially change the prediction made in pipe flow 
from simple 2 dimensional streamwise rolls to a spanwise and cross-stream localised 3 dimensional 
state. This new nonlinear optimal is demonstrably more efficient in triggering turbulence than the 
linear optimal indicating that there are better ways to design perturbations to achieve transition. 

PACS numbers: 



Shear flows are ubiquitous in nature and engineering, 
and understanding how and why they become turbulent 
has huge economic implications. This has led to a num- 
ber of simplified canonical problems being studied such 
as plane Couette flow, channel flow and pipe flow which 
commonly exhibit turbulent behavior even when the un- 
derlying laminar state is linearly stable. In this case, a 
finite amplitude perturbation is required in order to trig- 
ger turbulence and a leading question is then what is the 
'most dangerous' or 'smallest' such perturbation (with 
the metric typically being energy). Beyond its intrinsic 
interest, such information is fundamentally important for 
devising effective control strategies to delay the onset of 
turbulence. 

Linear transient growth analysis has commonly been 
used to suggest the structure of such dangerous distur- 
bances The basic premise being that disturbances 
which experience the most (transient) growth are also 
most efficient at modifying the underlying shear to pro- 
duce instability and probable transition to turbulence 
0, 3- Recent improvements in our understanding of 
the laminar-turbulent boundary or 'edge', which deter- 
mines whether a given initial condition leads to a turbu- 
lent episode or to the laminar state, has supported this 
idea albeit extended to disturbances of finite amplitude. 
This is because some regions of this boundary have a 
much smaller energy level (e.g. 3) than the attracting 
region of the boundary-confined dynamics [f| so that the 
minimum energy point on the edge (the most dangerous 
disturbance) must experience considerable energy growth 
as it sweeps up to the attracting plateau. 

However, there are two tacit assumptions in using lin- 
ear transient growth to identify critical disturbances: (1) 
that the energy growth experienced by the most dan- 
gerous disturbance is the largest (or near largest) pos- 
sible at the critical energy; and (2) the optimal distur- 
bance which emerges from linear transient growth analy- 
sis reasonably approximates the finite-amplitude optimal 
for the properly nonlinear growth calculation. Efforts to 
test these assumptions have concentrated on very small 



dimensional systems [6j or restricted the search for dan- 
gerous disturbances within small subspaces [3, 0| with 
broadly supportive results. The only study to specifically 
test assumption (2) used the Blasius approximation for 
the boundary layer Q and found no qualitative difference 
between the nonlinear and linear optimals. 

In this letter, however, we show for the first time us- 
ing the full Navier-Stokes equations how nonlinearity can 
fundamentally change the optimal which emerges from a 
transient growth analysis in pipe flow at subcritical en- 
ergy levels, thereby contradicting assumption (2). The 
significance of this new state is that: a) it provides a 
much more efficient way to trigger transition than the lin- 
ear optimal; and b) it is 3-dimensional and shows signs 
of localisation thereby appearing more physically rele- 
vant than the 2-dimensional streamwise-independent lin- 
ear optimal. 

The transient growth problem is the optimisation ques- 
tion: what initial condition u(x, t — 0) (added as a per- 
turbation to the laminar flow) for the governing Navier- 
Stokes equations with fixed (perturbation) kinetic energy 
Eq will give rise to the largest subsequent energy Et at 
a time t — T later. This corresponds to maximising the 
functional 
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where ( } represents volume integration; (s, 0, z) are 
cylindrical coordinates directed along the pipe; A, i/(x, t), 
LT(x, t) and T(t) are Lagrange multipliers imposing the 
constraints of initial energy Eq, that the Navier-Stokes 
equations hold over t G [0,T], incompressibility and con- 
stant mass flux in time respectively (the system has been 



2 



non-dimensionalised by the pipe diameter D and the bulk 
velocity U so that Re := pUD/p, where p is the density 
and p is the dynamic viscosity, and the laminar flow is 
2(1 — 4s 2 )z ). Vanishing of the variational derivatives re- 
quires that u must evolve according to the Navier-Stokcs 
equations, v evolves according to the adjoint- Navier- 
Stokes equations and at times t — and T we have 
optimality and compatibility conditions linking the two 
sets of variables (e.g. sec [8] for details of the linearised 
problem). The method of solution is one of iteration as 
follows: 

• Make an initial guess for u(x, t = 0) and allow the 
flow to evolve according to the Navier-Stokes equa- 
tions until t = T . 

• Solve the compatibility condition for i/(x, T), 
<LSf/<5u(x,T) ee u(x,T) - f(x,T) = 0. 

• Allow the incompressible field i/(x, t) to evolve 
backwards in time until t = via the adjoint- 
Navier-Stokes equations 
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• Move u(x, 0) in the direction of the variational 
derivative 53?/5u(x, 0) ee -Au(x, 0) + i>(x,0) to 
increase «5f and repeat. 

The algorithm should converge if E does not exceed the 
critical energy for transition. Beyond this, sensitivity 
to initial conditions when u(x, T) reaches the turbulent 
state will lead to non-smoothness. 

Both direct and adjoint equations were solved using 
a fully spectral, primitive variables approach. Time 
stepping was done using a second order fractional step 
scheme, checked carefully against the code of 9]. The 
computational domain was a short periodic domain of 
length 7r radii with typical spatial resolution of 29 real 
Fourier modes azimuthally, 11 real Fourier modes axially 
and 25 modified Chebyshev polynomials radially in each 
of the 8 scalar fields (it, v, w,p, z/i, 1/3,11). All results 
have been checked for robustness to resolution changes. 
Retention of the nonlinear terms poses a fresh technical 
challenge: although the adjoint equation is linear in u, 
it is dependent on the evolution history of the forward 
variable u which now must be stored. 

The linear transient growth optimal u;i„(x;i?e) in 
pipe flow is well-known to be streamwise-independent 
(2D) rolls which evolve into much larger streamwise- 
independent streaks [l[: see figs 1 and 2. Maximum 
growth occurs at T Un » 12.2 x i?e/1000 (D/U) [13]. In- 
troducing nonlinearity (ie increasing Eq from 0), setting 
T = 12.2 x Re/1000 (D/U) and allowing only 2D flows, 
leads smoothly to a modified 2D optimal U2£>(x; Eo,Re) 
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FIG. 1: The evolution of the linear and nonlinear optimals at 
Re — 1750. The blue (upper) line corresponds to the nonlin- 
ear optimal for Eq — 2 x 10~ J while the red (lower) line is the 
linear optimal (Eq — > 0). The nonlinear result produces more 
growth and actually reaches its maximum at a slightly earlier 
time than T. 



with monotonically decreasing growth consistent with 
previous simulations 2]. Opening the optimisation up 
to fully 3D flows initially just recovers the 2D result but 
once Eq crosses a small threshold E^d (1.35 x 10~ 5 at 
Re = 1750), a completely new optimal U3£>(x; Eq, Re) 
appears. This 3D optimal emerges from the optimisa- 
tion procedure after it initially appears to converge to 
the 2D optimal and then transiently visits an interme- 
diate state: see fig 3. Identifying this 'loss of stability' 
of the 2D optimal provided an efficient way to compute 
Ezr>(Re). All optimisation results were robust over three 
very different choices of starting flow: a) uu n with noise; 
b) the asymmetric travelling wave 11] ; and c) a turbu- 
lent flow snapshot (all rescaled to the appropriate initial 
energy). This supported our supposition that the algo- 
rithm samples all possible flows of a given energy to select 
the global optimiser although no proof is available. 

Given the intensity of the runs (O(200) iterations and 
each iteration requires integrating forwards and back- 
wards over the period [0,T]), one other Reynolds num- 
ber, Re = 2250, was selected to confirm our findings. 
Here, the new 3D state becomes the nonlinear optimal 
at E3D — 4.8 x 10~ 6 and has essentially the same ap- 
pearance as at Re — 1750: see fig. [2l Unlike the linear 
optimal which is globally simple in form and undergoes 
an evolution that is well established (rolls advecting the 
mean shear to generate streaks), the 3D optimal is lo- 
calised to one side of the pipe and initially has both rolls 
and streaks of comparable amplitude. Figures Q] and [2] 
show a new 2-stage evolution: a preliminary phase when 
the flow delocalises followed by a longer growth phase 
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FIG. 2: Three snapshots of the linear optimal (top) and five 
snapshots (middle & bottom) of the 3D optimal for Re — 1750 
and Eo = 2 x 10 during its evolution. Labels refer to figure 
[1] arrows indicate cross-sectional velocities and colours axial 
velocity beyond the laminar flow (white/light for positive and 
red/dark for negative: outside shade represents zero). The 
bar chart shows the ratio of energy in each streamwise Fourier 
mode of the initial nonlinear optimal (a) . 
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FIG. 3: Re = 1750, E = 2 x 10~ 5 . The iterations are seeded 
with a noisy version of the 2D optimal which converges to the 
3D optimal by way of an intermediate 'saddle' state (shown). 



where the flow structure stabilises to essentially two 2D 
large-scale slow streaks sandwiching one fast streak near 
the boundary. 

For E a > 2x 1CT 5 at 1750 and 6.25 x 1CT 6 at 2250, 
the iterative procedure fails to converge. In either case, a 
direct numerical simulation starting with the 3D optimal 
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FIG. 4: Re=2500. The green (lowest at t = 50 D/U) line 
shows the evolution of the 3D optimal when given initial en- 
ergy S c . Because it is on the laminar-turbulent boundary two 
nearly identical initial conditions diverge after, in this case, 
22QD/U. The blue (middle at t = 50 D/U) line is the evo- 
lution of the 2D optimal for the exact initial energy S l a m for 
which the streaks become linearly unstable. The red (upper 
at t = 50 D/U) line shows the 2D optimal given initial en- 
ergy Sl m and allowed to evolve until it reaches a maximum 
amplitude whereupon 0.1% by amplitude unstable perturba- 
tion is added. Again the laminar-turbulent boundary can be 
identified. 



at the highest energy value yielding convergence does not 
reveal a turbulent episode. This implies that the critical 
energy level, <§ c (Re), for transition has not been reached 
12J. The reasons for this energy 'gap' are unclear and 
leaves open the possibility that a further new optimal 
may emerge. It is worth noting that the end state of the 
3D nonlinear optimal decays more quickly that of the 2D 
optimal (which is in fact a close approximation of the 
least decaying eigenmode). Increasing the optimisation 
time, T, will therefore lead to the recovery of the 2D 
optimal at larger values of Eo . We then expect to able to 
converge onto the 3D optimal at correspondingly larger 
values of Eq, allowing the gap from E^d to S c to be 
closed. 

What can be tested, however, is whether the 3D op- 
timal is more efficient at triggering turbulence than the 
linear optimal when rescaled. Taking the initial condi- 
tion Au 3 fl(x;2 x 10~ 5 , 1750), we gradually increase the 
rescaling factor A until Eq = £ c {Re) is reached. Cal- 
culating the corresponding quantity for the linear opti- 
mal turns out to be less clearly denned because some 
3D noise is needed to trigger turbulence. As a result we 
make 2 different estimates, one strictly conservative and 
the other more realistic. The first <§ l s m is obtained by 
taking Auu n (x; 1750) and finding the initial energy for 
which the resultant streaks are just linearly unstable in 
this periodic domain @, Q. In the second <^ m , the same 
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FIG. 5: The effect of the initial energy on the growth of Auu n 
(red) and iu 3D (x,2 x 10 _5 ,1750) (blue) at Re = 1750. For 
small Eo the 2D result is the optimal but after E^d = 1-35 x 
10 -5 , the 3D optimal takes over. The vertical dashed line 
corresponds to S c , with the dotted lines being the relevant 
errorbars. Inset: The dependence of E$d (red), S c (blue) 
and g\ iTl (green) on Re (<fj in is even higher). For Re < 2000, 
errorbars on $ c indicate the energy range over which short to 
to extended turbulent episodes are triggered. 



approach (albeit with assumptions) for identifying criti- 
cal disturbances beyond exhaustive searches over initial 
conditions. As a result, the new 3D optimal found here 
supersedes the linear optimal as our current best theo- 
retical prediction for the most dangerous disturbance in 
pipe flow. 

There are two key directions for improving the result 
presented here: performing a further growth maximisa- 
tion over T and adopting larger, more realistic flow do- 
mains. Both represent formidable extensions even with 
today's computing power. After all, the discovery of the 
first true nonlinear optimal has had to wait almost two 
decades after the linear result was established in pipe 
flow. In larger domains, we expect further localisation of 
the nonlinear optimal since energy is defined as a global 
quantity whereas nonlinearity is important wherever the 
velocity field is locally large. This strongly suggests that 
in a long pipe the optimal should localise fully (i.e. in 
the axial direction as well) which would make it an in- 
teresting focus for experiments. 

We thank the referees for their comments. The calcu- 
lations in this paper were carried out at the Advanced 
Computing Research Centre, University of Bristol. 



initial condition was used but 0.1% of the most unstable 
perturbation (as found from the previous computation) 
is added to the streaks when they reach maximum am- 
plitude. S l ™ should be a (low) conservative estimate but 
even this is O(10) times larger than S c at Re = 2500 - 
see figure H]- whereas the more realistic §]? n is O(100) 
times larger. 

In figure [Sf inset) we plot E 3 d, <%c and S\ m as a 
function of Re which emphasizes that the 2D optimal 
(for which the linear result is an excellent approxima- 
tion) ceases to be a global maximum at an energy (at 
least) several orders of magnitude before it approaches 
the laminar-turbulent boundary. The 3D optimal, in 
contrast, crosses the laminar-turbulent boundary only 
shortly after it emerges at E^p (e.g. at k 5 x 10~ 5 where 
E 3D = 1.35 x 10~ 5 at Re = 1750). This means that the 
energy growth experienced by the 3D optimal must in- 
crease dramatically with E which is illustrated in figure 
[5] at Re = 1750: this growth is now a lower bound on the 
maximum possible for Eq > 2 x 10~ 5 . Assuming that 
a 3D optimal will always appear at subcritical energies 
(reasonable as the 2D disturbance cannot trigger turbu- 
lence), the critical energy must be bounded from below 
by Ezd and above by S c . 

In this letter we present the first demonstration that in- 
cluding nonlincarities in the problem of transient growth 
substantially changes the form of the optimal at ener- 
gies below that needed to trigger turbulence. The sig- 
nificance of this result comes from the fact that tran- 
sient growth analysis is currently the only constructive 
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